Energy transfer and dual cascade in kinetic magnetized plasma turbulence 
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The question of how nonhnear interactions redistribute the energy of fluctuations across available 
degrees of freedom is of fundamental importance in the study of turbulence and transport in mag- 
netized weakly collisional plasmas, ranging from space settings to fusion devices. In this letter, we 
present a theory for the dual cascade found in such plasmas, which predicts a range of new behavior 
that distinguishes this cascade from that of neutral fluid turbulence. These phenomena are explained 
in terms of the constrained nature of spectral transfer in nonlinear gyrokinetics. Accompanying this 
theory are the flrst observations of these phenomena, obtained via direct numerical simulations us- 
ing the gyrokinetic code AstroGK. The basic mechanisms that are found provide a framework for 
understanding the turbulent energy transfer that couples scales both locally and non-locally. 



It is known that a dual cascade occurs in 2D neutral 
fluid turbulence [T] and also in some reduced models of 
magnetized plasma turbulence such as the Hasegawa- 
Mima equation [2]. This cascade has the intriguing fea- 
ture that two invariants are spectrally transfered in op- 
posite senses; that is, the enstrophy flows to small scales 
while the energy is organized at large scales. This be- 
havior is due to the fact that the redistribution of these 
two invariants is constrained by the relationship between 
their spectral densities, as first explained by Fj0rtoft [3]; 
also see [H [5] for application relevant to plasmas. 

In gyrokinetics, the turbulent field is a scalar distribu- 
tion function g |12| . which varies over both position and 
velocity space, i.e. the particle phase space. There are 
two important quadratic quantities associated with this 
description, namely the free energy W and the electro- 
static energy E, which are "nonlinear invariants" (that 
is, conserved by the sole action of the nonlinearity) in 
both 2D and 3D. These quantities represent the energy of 
turbulent fluctuations. To study their transfer, we focus 
on a minimal form of gyrokinetics which retains nonlinear 
dynamics, but neglects linear drive and damping mecha- 
nisms: the homogeneous 2D gyrokinetic equation |6j. 

It was theoretically predicted [6l [7] and numerically 
observed [5] that 2D gyrokinetics has a cascade which 
exhibits the forward transfer of free energy to smaller 
scales in both position and velocity space, while electro- 
static energy is transfered inversely to larger scales. In 
this letter, we find that this conventional dual cascade is 
only one possible behavior. Adapting the arguments of 
Fj0rtoft, a full range of behavior can be understood by 
considering how the relationship between the nonlinear 
invariants constrain evolution of the phase-space spec- 
trum. This analysis reveals that the direction and local- 
ity of the nonlinear transfer actually depend on the way 
in which the free energy is initially distributed (or in- 
jected) across the phase space spectrum. For certain ini- 
tial conditions, we find that a dramatic non-local inverse 
transfer of electrostatic energy will occur, whereby large 
scales are driven directly by small scales. In the other ex- 
treme, we find that the inverse transfer is replaced by a 
forward transfer which can be nonlocal or local. Where 



the transfer is nonlocal, this provides a mechanism for 
the damping of large scales by small scales, which could 
be exploited to reduce the turbulent transport that limits 
the performance of fusion devices. 

In 2D, the gyrokinetic distribution function g depends 
on the gyro-center position R = xi?^; -I- yRy, the per- 
pendicular velocity v and the time t. This field evolves 
nonlinearly via the E x B velocity determined from the 
electrostatic potential (p which is a function of position r 
and time t {i.e., see Eqn. (2.5) of [3]). 

A crucial assumption at the outset is the existence of 
an external scale in velocity space, which we take to be 
the thermal velocity wth- This determines the characteris- 
tic spatial scale, the (thermal) Larmor radius p — Uth/f^- 
We assume that g{v) is attenuated above t;th- In fact, for 
the perturbed entropy to be a finite quantity, g{v) must 
fall off at least as fast as a gaussian as u ^- cx). 

We normalize velocity and spatial quantities to the 
thermal velocity and Larmor radius, respectively (see [6] 
for the full normalization). In what follows we will focus 
on interactions among "sub-Larmor" scales, i.e. spatial 
scales and velocity scales that are smaller than the Lar- 
mor radius and thermal velocity, respectively. This is 
distinct from the typical fiuid limit fc <C 1, where the 
formulation of the energetics depends on the reduced de- 
scription specific to the particular problem of interest. 

To obtain a discrete spectral representation, we assume 
the system to have finite extent in both position and 
velocity space. Adapting the continuous representation 
from [5], we decompose g using a Bessel- Fourier series: 
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where pj = Xj/Vcut and A^ is the jth zero of the zeroth 
order Bessel function Jq and 9 is the Heaviside step func- 
tion, which explicitly provides the velocity cutoff. The 
cutoff velocity is t'cut(^i) ^ Ij chosen for each ki such 
that Pj = ki for some j. Inverting the transform gives 
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where L is the system size in R^ and Ry. In spectral 
form, the quasi-neutrahty constraint, which relates the 
potential to the distribution function, is simply 



<^(ki) =/3(ki)5(ki,fcj) 



(3) 



where /3 — 27r/(l + r — fg), r is a constant related to 
the electron response [S] and Fq ~ ^o(fc^)e^'^ where /q 
is the zeroth-order modified Bcssel function. For fc ^ 1, 
the regime of interest, Fq oc 1/k may be neglected. The 
electrostatic energy is defined 
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where, from Eqn. [s) E{\ii) = 7r/3|5(kj, fc,)p. The other 
invariant, being a function of velocity space, we call "ki- 
netic free energy": G{v) — 1/(2L^) Jd^R g"^. For the 
purposes of this paper, we are not interested in this 
full invariant but only on the consequences its existence 
has for constraining the spectral evolution of fluctuations 
(note that [9] do retain the details of G). The relevant 
quantity is formed by integrating G over its velocity- 
space dependence; henceforth, we focus on this quantity, 
the (generalized) free energy: 
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where W{'ki,pj) = 27r^pj|g(k,j,pj)p/(l + t). For com- 
parison we note that W is proportional to the quantity 
Wgi from [3]. Note that W (or more generally G{v)) and 
E, being the known quadratic invariants, play a special 
role in constraining spectral transfer, just as enstrophy 
and energy do in 2D fluid turbulence |13) . 

We now arrive at the central equation of this letter, a 
constrains on the spectral densities of the invariants. In 
the limit fc ^ 1, it follows from the above definitions that 
(see also the derivation "W , Eqn. (7.14)) 
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Upon initial inspection, this relationship is similar to 
that for the NS equation: ZNs(ki) = fc^i<JNs(ki), where 
-Ens and Zns are the energy and enstrophy, respectively. 
However, this relationship constrains energetic interac- 
tions involving those components of free energy residing 
only along the diagonal in k-p space (fc^ = pj) [Tl|. This 
fact underscores the difference between the roles of ve- 
locity space and position space and, as we will see, is the 
central reason for the novel behavior that is discovered. 
Let's first consider what Eqn. [6] implies for three-scale 
interactions. 
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(a) I: Fj0rtoft-type 



A: 
(b) II: Kinetic-type 



FIG. 1: Constrained energetic transitions involving 
three scales. The dotted diagonals indicate k ^ p. 



We define the free energy at the wavenumber pair 
{ki,pi) to be Wi = Wiki,pi) = E|k,|=fci ^(ki,pO 
and the corresponding electrostatic energy to be Ei — 
5{ki - pi)E{ki) = d{ki - Pi) S|kj|=fc, -^(^j)' where 6 is 
the discrete delta function. Let's consider the difference 
in energies between the initial and final times, to and 
ti, for three components i = 1, 2 and 3: Thus we de- 
fine AW, = W^ih) - W.ito) and AE, = E,{h) - E,{to). 
Conservation of W and E implies: 



AWi + AW2 + AW3 = 0, 
AEi + AE2 + AE3 ^ 0. 



(7) 



As illustrated in Fig. [T] there are two types of con- 
strained transfer. For definiteness, we take ki < k2 < k^. 

In type I transitions. Fig. |la[ all three components are 
diagonal, pi — ki. In this case, transfers are constrained 
in precisely the sense of Fj0rtoft, with the substitution 
kf — > ki. From equations Eqns.^andJT^we obtain: 
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- A:i)/(fc3 - fci). 


(9) 



The quantities in parentheses are positive (because ki < 
^2 < ^3)- Thus, as noted by Fj0rtoft, it is only the 
intermediate wavenumber, fc2, which can be a source for 
both the two remaining components. 



Type II transitions. Fig. lb involve two diagonal com- 
ponents, and also lead to inverse transfer of E for the 
arrow directions indicated. To see this, first note that 
AE3 = because ps ^ k^. From Eqns.[6]and[7]we find 



AEi = AW3/{k2 ~ fci). 



(10) 



Thus, a transfer of free energy from diagonal to non- 
diagonal components must be accompanied by a simul- 
taneous inverse transfer of E; as we will shortly see, the 
reverse process can also spontaneously occur. 

Finally, note that transitions involving only one diag- 
onal component are forbidden as they cannot conserve 
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FIG. 2: Example spectral distributions 
W{k,p)/W] and initial evolution. Diagonals 
marked by dotted lines. 



E, "while those that involve no diagonal components are 
unconstrained. 

The elementary three-scale transitions form building 
blocks for thinking about more complicated evolution. 
For instance, the for"ward diffusion of free energy along 
and away from the diagonal, generally speaking, must 
induce an inverse transfer of electrostatic energy. On 
the other hand, we can imagine that the free energy 
could be initially concentrated in non-diagonal compo- 
nents and, in this case, the random redistribution of free 
energy would excite diagonal components causing elec- 
trostatic energy to flow forward, in the reverse sense of 
Fig. |lb[ From these considerations, it seems that the 
tendency for electrostatic energy to flow forward or back- 
ward should depend on the way in which the free energy 
is distributed relative to the diagonal. We will use the 
parameter k = W/E to help quantify this: small k corre- 
sponds to free energy being concentrated on the diagonal 
while large k corresponds to free energy being concen- 
trated off-diagonal. 

Let us now turn to the numerical results. We use the 
nonlinear gyrokinetic code AstroGK jlOj, including col- 
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FIG. 3: Evolution of the electrostatic energy spectrum. 



lisional dissipation as a free energy sink at large k to 
prevents nonphysical "bottle-necking". We also restrict 
attention to cases where the lowest wavenumber is unex- 
cited, so as to avoid effects of energetic condensation at 
large scales. 

Having explored a variety of initial spectral distribu- 
tions, we present three here, corresponding to different 
magnitudes of k, which exemplify different types of spec- 
tral evolution. To accompany the discussion, we provide 
some simple analytic calculations. Fig. [2] shows the ini- 
tial evolution of W{k,p) for the three cases, in the order 
of increasing k. 

The first case. Fig. 2a]|2b is an extreme case, corre- 
sponding to small K. The free energy is concentrated 
around a single diagonal component (fco, kg). The trans- 
fer of free energy away from this initial peak can only be 



accomplished by either Fj0rtoft-type (Fig. la) or kinetic- 
type transitions (Fig. |lb[) , which both demand an inverse 



transfer of electrostatic energy. As shown in Fig. 2b this 
is indeed what occurs: diagonal components are strongly 
excited at wavenumbers smaller than fco- As seen in 
Fig. [3aj the transfer that occurs is actually nonlocal, with 
a dominant peak of electrostatic energy E{k) emerging 
spontaneously at a wavenumber, k^, that is significantly 
smaller than kg. 



Notice in Fig. 2b that the W{k,p) spectrum is excited 
along lines in the k-p plane. This hints to the presence 
of an instability that involves the coupling of a small 
number of Fourier components. We confirm that a four- 
wave truncation (see, e.g., |llj ) is in fact an excellent 
approximation for the instability that is induced in the 
presence of a single Fourier component; the details of this 
calculation will be reported elsewhere [13]. The growth 



rate curve (7s (fc) in arbitrary units) of this instability 



is incorporated into Fig. 3a Note that the wavenumber 
of peak growth rate is close to that of the peak which 
develops in the energy spectrum. 

At a later stage, the case of Fig. [2a| evolves to a state in 
which the free energy is more broadly distributed about 
the diagonal. In this state, the 4-wave instability is ap- 
parently suppressed, or otherwise lost in the bath of fluc- 
tuations. Thus, many more degrees of freedom partici- 
pate and the evolution becomes much more cascade-like, 
whereby the free energy diffuses locally in k-p space. 

To treat the cases of Figs. 2cpf more quantitatively 
we perform a calculation, extending the arguments of 
Fj0rtoft with some simple assumptions. In each case 
that we consider, the initial spectrum is dominated by 
a single wavenumber fco- During a finite period of time, 
energy is liberated from this distribution in the amount 
Wo and Eo, satisfying Wq/Eq = k. Some amount of 
these energies, W and E' , diffuse forward in k-p space. 
We observe that this redistribution is locally diffusive in 
k-p space, spreading randomly to unoccupied wavenum- 
bers. By simple mode counting, this implies the relation- 
ship W' w Ck^E' , where C is some (possibly universal) 
constant of order unity. Now, in order to maintain the 
balance of both invariants, there must be some transfer of 
electrostatic energy E^ to a wavenumber fc^, along with 
the corresponding free energy Wd = k^Ed- By conserva- 
tion of energies we have Wq — W'-\-Wd and Eq = E'+Ed- 
We can combine these expressions to obtain 



R={n-kd)l{Ckl-K), 



(11) 



where R ~ E' / Ed- The parameters R and kd characterize 
the direction and locality of the transfer of E. From 
Eqn. 11 it is clear that k is an important parameter: at 



the critical case Kc = Cfcg, R diverges and changes sign. 

In Fig. [2cj we have an initial condition with free en- 
ergy distributed broadly around k^ = 20. This case is 
the closest in behavior to the conventional dual cascade 
of fluid mechanics. The free energy is observed to dif- 
fuse locally forward in k-p space whereas the transfer of 
electrostatic energy is mostly local and directed in the 
inverse direction, to k < ko, as is apparent in Fig. |3b| 
Thus it appears that < i? < 1, so, since k > k^ > kd, 
we can conclude that k < Kc- 

For K ^ Kc — C/cq, the evolution of E{k) changes 
character sharply. If, for instance, the magnitude of kd is 
comparable to fcg, i.e. kd ~ O{ko), then, by Eqn. 11 R 
is a large positive number for k '^ k^. Therefore most of 
the electrostatic energy transfered is forward. Now, if k 
exceeds Kc, R becomes negative and electrostatic energy 
is extracted from the kd component. For kd > kg, this 
then corresponds to an inverse transfer of E, while for 
kd < fco; the transfer is forward. 

In Figs. [2ep^ we have the case k ^ Kc- Here there 
is a dominant bath of free energy around A: ~ 20 while 



we fix kd at the outset by providing an initial "seed" 
energy at kd ^ 5. As the system evolves, electrostatic 
energy is drawn from the component at kd and so, by 
Eqn. [TT] Ed is negative. This process must proceed by 
kinetic transitions in the reverse sense of Fig.[Tbj whereby 
electrostatic energy is "forcefully" drawn from the large 
scale components to support the diffusion of free energy 
at small scales. 

The local diffusion of free energy is clearly observed 
from an initial bath at high-fc and \ow-p. Also, the evo- 



lution of E(k), shown in Fig. 3c shows forward trans- 
fer of electrostatic energy. For comparison. Fig. [3d| 
shows another large-K case, but with the seed energy at 
kd — ko. Here the system evolves by a dual forward 
cascade, whereby both W and E are transfered locally 
to finer scales. Note that the surprising "cascade rever- 
sal" of E is essentially due to the freedom afforded by the 
larger spectral space available to W. Note that the direc- 
tionality and locality of transfers in all of the above cases 
has been confirmed by a nonlinear transfer diagnostic; 
these observations will be reported elsewhere. 

It is true that in a general 3D gyrokinetic system, 
magnetic geometry and linear mode physics can formally 
break the invariance of E and W by introducing both 
damping and instability; this strongly affects the nature 
of the turbulence. However, we stress that the nonlin- 
ear transfer of fiuctuation energy must always occur in 
such a way as to conserve both E and W. This is due to 
the fact that the nonlinearity can always be written in 
terms of a Poisson bracket involving a stream function (in 
the electrostatic limit, the electrostatic potential) for the 
turbulent fiow; this is not possible in 3D fluid turbulence 
(or general plasma turbulence) and thus energy transfer 
is not constrained as it is in 2D. Because of these facts, 
we believe that the mechanisms revealed in this letter are 
a generic feature of kinetic magnetized plasmas. 
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